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Abstract 

The link between a particular class of growth processes and ran- 
dom matrices was established in the now famous 1999 article of Baik, 
Deift, and Johansson on the length of the longest increasing subse- 
quence of a random permutation |BDJ99j . During the past ten years, 
this connection has been worked out in detail and led to an improved 
understanding of the large scale properties of one-dimensional growth 
models. The reader will find a commented list of references at the 
end. Our objective is to provide an introduction highlighting random 
matrices. From the outset it should be emphasized that this connec- 
tion is fragile. Only certain aspects, and only for specific models, the 
growth process can be reexpressed in terms of partition functions also 
appearing in random matrix theory. 
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1 Growth models 



A growth model is a stochastic evolution for a height function h{x,t), x 
space, t time. For the one-dimensional models considered here, either x G M 
or a; G Z. We first define the TASEF0 (totally asymmetric simple exclusion 
process) with parallel updating, for which h,x,t G Z. An admissible height 
function, h, has to satisfy h{x + 1) — h{x) = ±1. Given h{x, t) and x^, a local 
minimum of h{x,t), one defines 



independently for all local minima, and h{x,t + 1) = h{x,t) otherwise, see 
Figure [1] (left). Note that if h{-,t) is admissible, so is h{-,t + 1). 

There are two limiting cases of interest. At the waiting time for an 
increase by 2 has the geometric distribution (1 — q)q"', n = 0,1, . . .. Taking 
the g — )■ 1 limit and setting the time unit to 1 — g one obtains the exponential 
distribution of mean 1. This is the time-continuous TASEP for which, count- 
ing from the moment of the first appearance, the heights at local minima are 
increased independently by 2 after an exponentially distributed waiting time. 
Thus h,x & Tj and t G M. 

The second case is the limit of rare events (Poisson points), where the 
unit of space and time is and one takes g — )■ 0. Then the lattice spacing 
and time become continuous. This limit, after a slightly different height 
representation (see Section [5] for more insights), results in the polynuclear 
growth model (PNG) for which /i G Z, x, t G M. An admissible height 
function is piecewise constant with jump size ±1, where an increase by 1 
is called an up-step and a decrease by 1 a down-step. The dynamics is 
constructed from a space-time Poisson process of intensity 2 of nucleation 
events. h{x,t) evolves deterministically through (1) up-steps move to the 
left with velocity —1, down-steps move to the right with velocity +1, (2) 
steps disappear upon coalescence, and (3) at points of the space-time Poisson 
process the height is increased by 1, thereby nucleating an adjacent pair of up- 
step and down-step. They then move symmetrically apart by the mechanism 
described under (1), see Figure [1] (right). 

The TASEP and PNG have to be supplemented by initial conditions 
and, possibly by boundary conditions. For the former one roughly divides 
between macroscopically flat and curved. For the TASEP examples would 

^Here we use the height function representation. The standard particle representa- 
tion consists in placing a particle at x if h{x + 1) — h{x) = — 1 and leaving empty if 
h{x + 1) - h{x) = 1. 



h{x^,t + 1) 



h{x^:,t) + 2 with probability 1 
h{x^:,t) with probability g 



g, < g < 1 



(1) 
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Figure 1: Growth of TASEP interface (left) and continuous time PNG (right) 



be h{x) = for x even and h{x) = 1 for x odd, which has slope zero, 
and h{x + 1) — h{x) independent Bernoulli random variables with mean m, 
|m| < 1, which has slope m. An example for a curved initial condition is 
h{0) = 0, h{x + 1) - h{x) = -1 for x = -1, -2, . . ., and h{x + 1) - h{x) = 1 
for a; = 0, 1, . . .. 

What are the quantities of interest? The most basic one is the macro- 
scopic shape, which corresponds to a law of large numbers for 

lh{[yt],[st]) (2) 

in the limit t — )■ oo with ?/ G R, s G M, and [•] denoting integer part. From 
a statistical mechanics point of view the shape fluctuations are of prime 
concern. For example, in the flat case the surface stays macroscopically flat 
and advances with constant velocity v. One then would like to understand 
the large scale limit of the fluctuations {h{x, t) — vt, {x, t) G Z x Z}. As will 
be discussed below, the excitement is triggered through non-classical scaling 
exponents and non-Gaussian limits. 

In 1986 in a seminal paper Kardar, Parisi, and Zhang (KPZ) proposed 
the stochastic evolution equation |KPZ86j 

d f d \2 9^ 

—h{x,t) = \{^—h{x,t)j +uo — h{x,t)+T]{x,t) (3) 

for which h,x,t G M. f]{x,t) is space-time white noise which models the 
deposition mechanism in a moving frame of reference. The nonlinearity re- 
flects the slope dependent growth velocity and the Laplacian with uq > 
is a smoothing mechanism. To make the equation well-deflned one has 
to introduce either a suitable spatial discretization or a noise covariance 
{ri{x,t)r]{x' ,t')) = g{x — x')6(t — t') with g{x) = g{—x) and supported close 
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to 0. KPZ argue that, according to with initial conditions h{x, 0) = 0, 
the surface width increases as t^^^, while lateral correlations increase as 
It is only through the connection to random matrix theory that universal 
probability distributions and scaling functions have become accessible. 

Following TASEP, PNG, and KPZ as guiding examples it is easy to con- 
struct variations. For example, for the TASEP one could introduce evap- 
oration through which heights at local maxima are decreased by 2. The 
deposition could be made to depend on neighboring heights. Also general- 
izations to higher dimensions, x G M°' or x G Z'^, are easily accomplished. 
For the KPZ equation the nonlinearity then reads X(\/xh{x,t)^ and the 
smoothening is i>oAh{x,t). For the PNG model in d = 2, at a nucleation 
event one generates on the existing layer a new layer of height one consisting 
of a disk expanding at unit speed. None of these models seem to be directly 
connected to random matrices. 

2 How do random matrices appear? 

Let us consider the PNG model with the initial condition h{x, 0) = 
under the constraint that there are no nucleations outside the interval 
[— t > 0, which also refered to as PNG droplet, since the typical 
shape for large times is a semicircle. We study the probability distribu- 
tion F{{h{0,t) < n}), which depends only on the nucleation events in the 
quadrant {{x, s)\ \x\ < s, s — t < x < t — s , < s <t}. Let us denote by 
CO = {uj^^\ . . . , co^^'^) a set of nucleation events and by h{0,t]uj) the corre- 
sponding height. The order of the coordinates of the co^^^^s in the frame 
{x = ±t} naturally defines a permutation of n elements. It can be seen 
that h{0, t; co) is simply the length of the longest increasing subsequence of 
that permutation. By the Poisson statistics of u, the permutations are ran- 
dom and their length is Poisson distributed. By this reasoning, somewhat 
unexpectedly, one finds that F{{h{0,t) < n}) can be written matrix 
integral |BR01aj . Let Un be the set of all unitaries on C" and dU be the 
corresponding Haar measure. Then 



(jlj) can also be expressed as Fredholm determinant on £2 = hC^)- On £2 
we define the linear operator B through 



P({/i(0,t) < n 



}) = f 

Ju, 



d[/exp[ttr(t/ + f/*)]. 



(4) 



{Bf){x) 



/(x + 1) - /(x - 1) + -/(x) 



(5) 
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and denote by P<o the spectral projection onto B <0. Setting 9n{x) = 1 for 
X > n and 6n{x) = for x < n, one has 

P({/i(0, t) < n}) = det(l - e^P^oOn) . (6) 

Such an expression is famihar from GUE random matrices. Let Ai < . . . < Aat 
be the eigenvalues of an x GUE distributed random matrix. Then 

P({Aiv < y}) = det(l - OyPNOy) . (7) 

Now the determinant is over L^(]R, dy). If one sets H = —^-^ + jj^y'^, then 

P/v projects onto H < \/N. 

For large one has the asymptotics 

\j^^2N + N'/^^2 (8) 

with ^2 a Tracy- Widom distributed random variable, that is, 

< s}) = F^is) := det(l - XsKmXs), (9) 

with det the Fredholm determinant on L^(M, dx), Xsi^) = l(x > s), and K^i 
is the Airy kernel (see f[T^ below). Hence it is not so surprising that for the 
height of the PNG model one obtains [PSOObj 

h{0,t) = 2t + t'^^^2 (10) 

in the limit t — t- oo. In particular, the surface width increases as t^^^ in accor- 
dance with the KPZ prediction. The law in (ITOj) is expected to be universal 
and to hold whenever the macroscopic profile at the reference point, a; = in 
our example, is curved. Indeed for the PNG droplet h(x,t) = 1t^\ — {x/tY 
for large t, |x| < t. 

One may wonder whether ffTOl) should be regarded as an accident or 
whether there is a deeper reason. In the latter case, further height statistics 
might also be representable through matrix integrals. 

3 Multi-matrix models and line ensembles 

For curved initial data the link to random matrix theory can be understood 
from underlying line ensembles. They differ from case to case but have the 
property that the top line has the same statistics in the scaling limit. We first 
turn to random matrices by introducing matrix-valued diffusion processes. 

Let B{t) be GUE Brownian motion, to say B{t) is an A^ x A^ hermitian 
matrix such that, for every / G C^, t i— )■ {f,B(t)f) is standard Brownian 
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motion with variance (/, /)^ min(t, t') and such that for every unitary U it 
holds 

UB{t)U* = B{t) (11) 

in distribution. The N x N matrix- valued diffusion, A{t), is defined through 
the stochastic differential equation 

dA{t) = -V'{A{t))dt + dB{t) , A{0) = A, (12) 

with potential \/ : M M. We assume A = A*, then also A{t) = A{t)* for all 
t > 0. It can be shown that eigenvalues of A{t) do not cross with probability 
1 and we order them as Ai(t) < ... < XN{t). t ^ (Ai(t), . . . , Aiv(i)) is the 
line ensemble associated to f lT^ . 

In our context the largest eigenvalue, A7v(^) is of most interest. For 
— )■ oo its statistics is expected to be independent of the choice of V . We 
first define the hmit process, the Airy2 process ^2(^)5 through its finite di- 
mensional distributions. Let 

^Ai = -^ + 2/ (13) 

as a self-adjoint operator on L^(R, dy). The Airy operator has M 
as spectrum with the Airy function Ai as generalized eigenfunctions, 
HxiPd{y — A) = AAi(y — A). In particular the projection onto {H^i < 0} 
is given by the Airy kernel 

KAiiy, y') = / dAAi(|/ + A)Ai(|/' + A) . (14) 

The associated extended integral kernel is defined through 

KA,iy,T;y',T') = -ie-^^'--^"--'){y,y')liT' > r) + (e^^-i^Aie-^'^^Od/, 1/0- 

(15) 

Then, the m-th marginal for A2{t) at times ti < t2 < • • • < is expressed 
as a determinant on L^(]R x {ti, . . . ,tm}) according to |PS02bj 

^{A2{tl) < Si,...,A2{tm) < Sm) = det{l - XsKa.Xs) , (16) 

with Xs{x,ti) = l{x > Si), t I— )■ A2{t) is a stationary process with continuous 

sample paths and covariance (72 (^) = Cov(^2(0), ^2(^)) = Var(^2) — |i^|+C^(^^) 
for t -> and g2{t) = f-^ + C(r^) for |t| oo. 

The scaling limit for A2{t) can be most easily constructed by two slightly 
different procedures. The first one starts from the stationary Ornstein- 
Uhlenbeck process A^^{t) in ([I2D, which has the potential V{x) = x^/2N. 
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Its distribution at a single time is GUE, to say Z]^^ exp[—j^tTA'^]dA, where 
the factor results from the condition that the eigenvalue density in the 
bulk is of order 1. Let A°^(t) be the largest eigenvalue of A^^(t). Then 

lim Ar-i/3(AOU(2iV2/3t) - 2N) = A2{t) (17) 

in the sense of convergence of finite dimensional distributions. The N'^/^ 
scaling means that locally A^^(t) looks like Brownian motion. On the other 
side the global behavior is confined. 

The marginal of the stationary Ornstein-Uhlenbeck process for two time 
instants is the familiar 2-matrix model |EM98l[NF98] . Setting = A°^{0), 
A2 = A'~'^(t), t > 0, the joint distribution is given by 

^exp ( - ^^^^^L__tT[Al + Al- 2gAiA2])dAidA2, q = exp{-t/2N). 

^ (18) 
A somewhat different construction uses the Brownian bridge defined by 
([T^ with y = and A^^{-T) = A^^{T) = 0, that is 

A^^it) = B{T + t)-'^B{2T), \t\<T. (19) 

The eigenvalues t 1— )■ (A^^(t), . . . , A^^(t)) is the Brownian bridge line ensem- 
ble. Its largest eigenvalue, X^{t), has the scaling limit, for T = 2N, 

lim (Xl^{2N^/H) - 2N) +f = A(t) (20) 

in the sense of finite-dimensional distributions. Note that A^^(t) is curved 
on the macroscopic scale resulting in the displacement by —t^. But with this 
subtraction the limit is stationary. 

To prove the limits f lT7|) and fl20|l one uses in a central way that the 
underlying line ensemble are determinantal. For the Brownian bridge this can 
be seen by the following construction. We consider iV independent standard 
Brownian bridges over [— T, T], &^^(t), j = 1,...,A^, 6j^(±T) = and 
condition them on non-crossing for |t| < T. The resulting line ensemble 
has the same statistics as A™(t), \t\ < T, j = 1,...,N, which hence is 
determinantal. 

The TASEP, and its limits, have also an underlying line ensemble, which 
qualitatively resemble {A^^(t), j = 1, . . . , N}. The construction of the line 
ensemble is not difficult, but somewhat hidden. Because of lack of space 
we explain only the line ensemble for the PNG droplet. As before t is the 
growth time and x is space which takes the role of t from above. The top 
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line is Xo{x,t) = h{x,t), h the PNG droplet of Section [T] Initially we add 
the extra lines Xj{x, 0) = j, j = —1, —2, .... The motion of these lines is 
completely determined by h{x,t) through the following simple rules: (1) and 
(2) from above are in force for all lines Xj{x, t), j = 0, —1, .... (3) holds only 
for Xo{x, t) = h{x, t). (4) The annihilation of a pair of an adjacent down-step 
and up-step at line j is carried out and copied instantaneously as a nucleation 
event to line j — 1. 

We let the dynamics run up to time t. The line ensemble at time t 
is {Xj{x,t),\x\ < t,j = 0,-1,...}. Note that Xj{±t,t) = j. Also, for a 
given realization of h{x,t), there is an index jo such that for j < jo it holds 
Xj{x, t) = j for all X, \x\ < t. The crucial point of the construction is to have 
the statistics of the line ensemble determinantal, a property shared by the 
Brownian bridge over [— t, t], Xj^{x), \x\ < t. The multi-line PNG droplet 
allows for a construction similar to the Brownian bridge line ensemble. We 
consider a family of independent, rate 1, time-continuous, symmetric nearest 
neighbor random walks on Z, rj{x), j = 0, —1, . . .. The j-th random walk 
is conditioned on rj{±t) = j and the whole family is conditioned on non- 
crossing. The resulting line ensemble has the same statistics as the PNG line 
ensemble {Xj{x,t),j = 0, —1, . . .}, which hence is determinantal. 

In the scaling limit for x = 0{t'^^^) the top line Ao(a;, t) is displaced by 2t 
and t^/^ away from X-i{x,t). Similarly A^^(x) for x = 0{N'^^^) is displaced 
by 2N and order N'^^'^ apart from Xf^^{x). On this scale the difference 
between random walk and Brownian motion disappears but the non-crossing 
constraint persists. Thus it is no longer a surprise that 

hmt-^/^(h(f^^x,t)-2t)+x^ =A2(x), xeR, (21) 

in the sense of finite dimensional distributions |PS02b] . 

To summarize: for curved initial data the spatial statistics for large t is 
identical to the family of largest eigenvalues in a GUE multi-matrix model. 

4 Flat initial conditions 

Given the unexpected connection between the PNG model and GUE multi- 
matrices, a natural question is whether such a correspondence holds also for 
other symmetry classes of random matrices. The answer is affirmative, but 
with unexpected twists. Consider the flat PNG model with h{x, 0) = and 
nucleation events in the whole upper half plane {{x,t),x G M, t > 0}. The 
removal of the spatial restriction of nucleation events leads to the problem 
of the longest increasing subsequence of a random permutation with involu- 
tion |BR01b(,[PS00bj . The limit shape will be flat (straight) and in the limit 
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t — )■ cxD one obtains 

/i(0,t) = 2t + el^'/^ (22) 

where 

P(6 < s) = Fi(2-2/3s). (23) 

The distribution function Fi is the Tracy- Widom distribution for the largest 
GOE eigenvalue. 

As before, we can construct the line ensemble and ask if the link to 
Brownian motion on GOE matrices persists. Firstly, let us compare the line 
ensembles at fixed position for flat PNG and at fixed time for GOE Brownian 
motions. In the large time (resp. matrix dimension) limit, the edges of these 
point processes converge to the same object: a Pfaffian point process with 
2x2 kernel |Fer04] . It seems then plausible to conjecture that also the 
two line ensembles have the same scaling limit, i.e., the surface process for 
fiat PNG and for the largest eigenvalue of GOE Brownian motions should 
coincide. Since the covariance for the fiat PNG has been computed exactly 
in the scaling limit, one can compare with simulation results from GOE 
multi-matrices. The evidence strongly disfavors the conjecture |BFP08] . 

The process describing the largest eigenvalue of GOE multi-matrices 
is still unknown, while the limit process of the fiat PNG interface is 
known |BFS08j and called the Airyi process, ^i, 

hm r^'^'ihie'^x.t) - 2t) = 21/3^1 (2-2/3a;). (24) 

Its m-point distribution is given in terms of a Fredholm determinant of the 
following kernel. Let B[y,y') = Ai{y + y'), Hi = — Then, 

KAAy,r;y',T') = -(e~(-'--)^^)(y,y')l(r' > r) + (e-^^i?e-^'^^)(y, y') (25) 

and, as for the Airy2 process, the m-th marginal for Ai{t) at times 
ti < t2 < • • • < is expressed through a determinant on L^(]Rx {ti, . . . , tm}) 
according to 

F(A(ti) < Si,...,A(U < Sm) =det{l -XsKa.Xs), (26) 

with Xs{x,ti) = l{x > Si) |Sas05llBFPS07j . The Airyi process is a stationary 
process with covariance gi(t) = Cov(^i(0), ^i(t)) = Var(^i) — \t\ + C(t^) for 
t — i- and gi{t) — t- super-exponentially fast as \t\ — )■ 00. 

The Airyi process is obtained using an approach different from the PNG 
line ensemble, but still with an algebraic structure encountered also in ran- 
dom matrices (in the the GUE- minor process |JN06|[OR06j ). We explain the 
mathematical structure using the continuous time TASEP as model, since 
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the formulas are the simplest. For a while we use the standard TASEP rep- 
resentation in terms of particles. One starts with a formula by Schiitz |Sch97j 
for the transition probability of the TASEP particles from generic positions. 
Consider the system of particles with positions xi(t) > X2(t) > ... > XN(t) 
and let 

G{xi, . . . , xn', t) = "^{xiit) 
Then 

G{xi, ...,XN;t) 

with 

Fnix,t) 

The function F„ satisfies the relation 

Fnix,t) = J2Fn-iiy,t). (30) 

y>x 

The key observation is that (125]) can be written as the marginal of a de- 
terminantal line ensemble, i.e. of a measure which is a product of determi- 
nants |Sas05] . The "lines" are denoted by with time index n, 1 < n < N, 
and space index i, 1 < i < n. We set Xi = Xn- Then 

Af-l 

G{xu ...,x^-t)= Yl ( n det{Mx7,x]^'))h=?l det{K-^{xf))Z=l 

x",2<i<n<N n=l 

(31) 

with ^'j^_i(x) = F_i+i{x-yN+i-i,t), 0„(x,?/) = l{x > y) and0„«+i,?/) = 1 
(here plays the role of a virtual variable). The line ensembles for the 
PNG droplet and GUE- valued Brownian motion have the same determinan- 
tal structure. However in (13T|) the determinants are of increasing sizes which 
requires to introduce the virtual variables x^^^. However, from the algebraic 
perspective the two cases can be treated in a similar way. As a result, the 
measure in fl3ip is determinantal (for instance for any fixed initial conditions, 
but not only) and has a defining kernel dependent on yi, . . . ,yiy. The dis- 
tribution of the positions of TASEP particles are then given by a Fredholm 
determinant of the kernel. To have flat initial conditions, one sets yi = N—2i, 
takes first the N ^ oo limit, and then analyzes the system in the large time 
limit to get the Airyi process defined above. 

A couple of remarks: 
(1) For general initial conditions (for instance for flat initial conditions), the 



xi, . . .,XN(t) = xn\xi{0) =yi,.. .,xn{0) = yn)- 

(27) 



det {Fi_j{xN+i-i - yN+i-j, t))i<i 



(28) 



(£) dw 



2Trijr^l>i w^-"+i(w - 1)" ' 



(29) 
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measure on {x^} is not positive, but some projections, like on the TASEP 
particles {x^}, defines a probability measure. 

(2) The method can be applied also to step initial conditions (xfc(O) = —k, 
k = 1,2, . . .) and one obtains the Airy2 process. In this case, the measure on 
{x"} is a probability measure. 

(3) In random matrices a measure which is the product of determinants of 
increasing size occurs too, for instance in the GUE-minor process |JN06l 



5 Growth models and last passage percola- 
tion 

For the KPZ equation we carry out the Cole-Hopf transformation 
Z{x,t) = exp{—Xh{x,t)/i'Q) with the result 

g^Zix,t) = - i-uo— + -r]ix,t)\ Z{x,t), (32) 

which is a diffusion equation with random potential. Using the Feynman-Kac 
formula, ( 132|) corresponds to Brownian motion paths, Xt, with weight 



(33) 



In physics this problem is known as a directed polymer in a random potential, 
while in probability theory one uses directed first/last passage percolation. 
The spirit of the somewhat formal expression (133|) persists for discrete growth 
processes. For example, for the PNG droplet we fix a realization, u, of the 
nucleation events, which then determines the height h{0,t;u) according to 
the PNG rules. We now draw a piecewise linear path with local slope m, 
\m\ < 1, which starts at (0,0), ends at (0,t), and changes direction only 
at the points of u. Let L{t; u) be the maximal number of Poisson points 
collected when varying over allowed paths. Then h{0,t;uj) = L(t;uj). So to 
speak, the random potential from fl55]) is spiked at the Poisson points. 

In this section we explain the connection between growth models and 
(directed) last passage percolation. For simplicity, we first consider the case 
leading to discrete time PNG droplet, although directed percolation can be 
defined for general passage time distributions. Other geometries like fiat 
growth are discussed later. 

Let oj{i,i), i,j > 1, be independent random variables with geometric 
distribution of parameter q, i.e., F{uj{i,j) = k) = {1 — q)q^, k > 0. An up- 
right path TT from (1, 1) to {n,m) is a sequence of points {it, with 
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Figure 2: Directed percolation and growth models 



- (ieje) e {(1, 0), (0, 1)}. The last passage time from (1,1) to 
{n, m) is defined by 



The connection between directed percolation and different growth models is 
obtained by appropriate projections of the three-dimensional graph of G. Let 
us see how this works. 

Let time t be defined hjt = i + j — 1 and position hj x = i — j. Then, 
the connection between the height function of the discrete time PNG and 
the last passage time G is simply |Joh03] 



Thus, the discrete PNG droplet is nothing else than the time-slicing along 
the i + j = t directions, see Figure [2J 

We can however use also a different slicing, at constant r = G, to obtain 
the TASEP at time r with step initial conditions. For simplicity, consider 
u}{i,j) to be exponentially distributed with mean one. Then u{n,m) is the 
waiting time for particle m to do his nth jump. Hence, G{n,m) is the time 
when particle m arrives at — m + n, i.e.. 



From the point of view of the TASEP, there is another interesting cut, 
namely at fixed j = n. This corresponds to look at the evolution of the 
position of a given (tagged) particle, x„(r). 




(34) 



h{x,t) = G{i,j). 



(35) 



F{G{n, m) <t) = P(a;„(r) >-m + n). 



(36) 
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A few observations: 

(1) Geometrically distributed random passage times correspond to discrete 
time TASEP with sequential update. 

(2) The discrete time TASEP with parallel update is obtained by replacing 
^ihJ) by + 1. 

The link between last passage percolation and growth holds also for gen- 
eral initial conditions. In f l34p the optimization problem is called point-to- 
point, since both (1, 1) and (n, m) are fixed. We can generalize the model by 
allowing uj{i,j) to be defined on G and not only for i,j > 1. Con- 
sider the line L = {i + j = 2} and the following point-to-line maximization 
problem: 

GL{n,m) = max N uj{i,j). (37) 

7r;L— >(n,m) ' ^ 
(i,j)eiT 

Then, the relation to the discrete time PNG droplet, namely 
h{x,t) = Giiiyj) still holds but this time h is the height obtained from flat 
initial conditions. For the TASEP, it means to have at time r = the parti- 
cles occupying 2Z. Also random initial conditions fit in the picture, this time 
one has to optimize over end-points which are located on a random line. 

In the appropriate scaling limit, for large time/particle number one ob- 
tains the Airy2 (resp. the Airyi) process for all these cases. One might wonder 
why the process seems not to depend on the chosen cut. In fact, this is not 
completely true. Indeed, consider for instance the PNG droplet and ask the 
question of joint correlations of h{x, t) in space-time. We have seen that for 
large t the correlation length scales as t^^^. However, along the rays leaving 
from {x,t) = (0,0) the height function decorrelates on a much longer time 
scale, of order t. These slow decorrelation directions depend on the initial 
conditions. For instance, for fiat PNG they are the parallel to the time axis. 
More generally, they coincide with the characteristics of the macroscopic sur- 
face evolution. Consequently, except along the slow directions, on the t^/^ 
scale one will always see the Airy processes. 

6 Growth models and random tiling 

In the previous section we explained how different growth models (TASEP 
and PNG) and observables (TASEP at fixed time or tagged particle motion) 
can be viewed as different projections of a single three-dimensional object. 
A similar unifying approach exists also for some growth models and random 
tiling problems: there is a 2 -f 1 dimensional surface whose projection to one 
less dimension in space (resp. time) leads to growth in 1-1-1 dimensions (resp. 
random tiling in 2 dimensions) |BF08a] . To explain the idea, we consider the 
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Figure 3: Illustration of the particle system and its corresponding lozenge 
tiling. The initial condition is on the left. 

dynamical model connected to the continuous time TASEP with step initial 
conditions, being one of the simplest to define. 

In Section H] we encountered a measure on set of variables 
Sn = {x"^, I < i < n < N}, see fl3T|) . The product of determinants of the 
(pnS implies that the measure is non-zero, if the variables belong to the 
set S*]^* defined through an interlacing condition, 

S'f = {x: G I x^' <x:< X-:,'}. (38) 

Moreover, for TASEP with step initial conditions, the measure on S'jy* is a 
probability measure, so that we can interpret the variable xf as the po- 
sition of the particle indexed by {i,n). Also, the step initial condition, 
Xn(t = 0) = 2;"(0) = —n, implies that xf{0) = i — n — 1, see Figure [3] (left). 

Then, the dynamics of the TASEP induces a dynamics on the particles 
in 5"]^* as follows. Particles independently try to jump on their right with 
rate one, but they might be blocked or pushed by others. When particle x^ 
attempts to jump: 

(1) it jumps if x" < x"^^, otherwise it is blocked (the usual TASEP dynamics 
between particles with same lower index), 

(2) if it jumps, it pushes by one all other particles with index {i + i,n + i), 
i > 1, which are at the same position (so to remain in the set 5]^*). 

For example, consider the particles of Figure [3] (right). Then, if in this state 
of the system particle (1, 3) tries to jump, it is blocked by the particle (1, 2), 
while if particle (2, 2) jumps, then also (3, 3) and (4, 4) will move by one unit 
at the same time. 

It is clear that the projection of the particle system onto {x", 1 < n < A^} 
reduces to the TASEP dynamics in continuous time and step initial condi- 
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t=2 t=l t=0 t=l t=2 




Figure 4: Two examples of configurations at time t = 2 obtained by the 
particle dynamics and its associated Aztec diamonds. The line ensembles 
are also drawn. 



tions, this projection being the sum in f lSTl) . 

The system of particles can also be represented as a tiling using three 
different lozenges as indicated in Figured The initial condition corresponds 
to a perfectly ordered tiling and the dynamics on particles reflects a corre- 
sponding dynamics of the random tiling. Thus, the projection of the model 
to fixed time reduces to a random tiling problem. In the same spirit, the 
shuffling algorithm of the Aztec diamond falls into place. This time the in- 
terlacing is ^Aztec = {^r I — ^i" — ^i+i} the dynamics is on discrete 
time, with particles with index n not allowed to move before time t = n. As 
before, particle (i, n) can be blocked by (i, n — 1). The pushing occurs in the 
following way: particle {i, n) is forced to move if it stays at the same position 
of (i — 1, 77, — 1), i.e., if it would generate a violation due to the possible jump 
of particle (i — l,n — 1). Then at time t all particles which are not blocked 
or forced to move, jump independently with probability q. As explained in 
detail in |Nor08] this particle dynamics is equivalent to the shuffling algo- 
rithm of the Aztec diamond (take g = 1/2 for uniform weight). In Figure S] 
we illustrate the rules with a small size example of two steps. There we also 
draw a set of lines, which come from a non-intersecting line ensemble simi- 
lar to the ones of the PNG model and the matrix- valued Brownian motions 
described in Section |3l 

On the other hand, let := — n. Then, the dynamics of the shuffling 
algorithm projected onto particles x", n > 1, is nothing else than the discrete 
time TASEP with parallel update and step initial condition! Once again, we 
have a 1 + 1 dimensional growth and a 2 dimensional tiling model which are 
different projections of the same 2 + 1-dimensional dynamical object. 
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7 A guide to the literature 



There is a substantial body of literature and only a rather partial list is 
given here. The guideline is ordered according to physical model under study. 

PNG model. A wide variety of growth processes, including PNG, are 
explained in |Mea98j . The direct link to the largest increasing subsequences 
of a random permutation and to the unitary matrix integral of [PS90] is 
noted in |PS00b|lPS00a] . The convergence to the Airy2 process is worked 
out in |PS02b] and the stationary case in |BR00|lPS04] . For flat initial con- 
ditions, the height at a single space-point is studied in |BR01a[lBR01b] and 
for the ensemble of top lines in |Fer04j . The extension to many space-points 
is accomplished by |BFS08j . Determinantal space-time processes for the 
discrete time PNG model are discussed by |Joh03] . External source at the 
origin is studied in |IS04b] for the full line and in |IS04a] for the half-line. 

Asymmetric simple exclusion (ASEP). As a physical model reversible 
lattice gases, in particular the simple exclusion process, were introduced 
by Kawasaki |Kaw72j and its asymmetric version by Spitzer |Spi70| . We 
refer to Liggett |Lig99| for a comprehensive coverage from the probabilistic 
side. The hydrodynamic limit is treated in |Spo91| and jKL99] . for example. 
There is a very substantial body on large deviations with Derrida as one 
of the driving forces, see [Der07j for a recent review. For the TASEP 
Schiitz |Sch97j discovered a determinantal-like formula for the transition 
probability for a finite number of particles on Z. TASEP step initial 
conditions are studied in the seminal paper by Johansson |JohOO] . The 
random matrix representation of [JohOOj may also be obtained by the Schiitz 
formula [NS04l[RS05j . The convergence to the Airy2 process in a discrete 
time setting |Joh05] . General step initial conditions are covered by |PS02aj 
and the extended process in |BFP09] . In |FS06bj the scaling limit of the 
stationary two-point function is proved. Periodic intial conditions were 
first studied by Sasamoto |Sas05j and widely extended in |BFPS07] . A 
further approach comes from the Bethe ansatz which is used in |GS92j to 
determine the spectral gap of the generator. In a spectacular analytic tour 
de force Tracy and Widom develop the Bethe ansatz for the transition proba- 
bihty and thereby extend the Johansson result to the PASEP [TWn9yTWn8] . 

2D tiling (statics). The connection between growth and tiling was first 
understood in the context of the Aztec diamond |JPS98j . who prove the 
arctic circle theorem. Because of the specific boundary conditions imposed, 
for typical tilings there is an ordered zone spatially coexisting with a 
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disordered zone. In the disordered zone the large scale statistics is expected 
to be governed by a massless Gaussian field, while the line separating the 
two zones has the Airy2 process statistics. 

a) Aztec diamond. The zone boundary is studied by |JPS98] and by |JN06] . 
Local dimer statistics are investigated in |CEP96] . Refined details are 
the large scale Gaussian statistics in the disordered zone |KenOO] . the 
edge statistics |Joh05] . and the statistics close to a point touching the 
boundary |JN06j . 

b) Ising corner. The Ising corner corresponds to a lozenge tiling under the 
constraint of a fixed volume below the thereby defined surface. The largest 
scale information is the macroscopic shape and large deviations |CK01] . 
The determinantal structure is noted in |UR03j . This can be used to study 
the edge statistics |FS03[[FPS04j . More general boundary conditions (skew 
plane partitions) leads to a wide variety of macroscopic shapes |UR07] . 

c) Six-vertex model with domain wall boundary conditions, as introduced 
in |KZJOO] . The free energy including prefactors is available jBF06] . A 
numerical study can be found in |AR05] . The mapping to the Aztec 
diamond, on the free Fermion line, is explained in |FS06aj . 

d) Kasteleyn domino tilings. Kasteleyn |Kas63j noted that Pfaffian methods 
work for a general class of lattices. Macroscopic shapes are obtained 
by |KQS06] with surprising connections to algebraic geometry. Gaussian 



fluctuations are proved in [KenOSj . 

2D tiling (dynamics), see Section El The shuffling algorithm for the Aztec 
diamond is studied in |EKLP92"|IPro03ilNor08j . The pushASEP is introduced 
by |BF08b] and anisotropic growth models are investigated in |BF08a] . see 
also |PS97] for the Gates- West cott model. A similar intertwining structure 
appears for Dyson's Brownian motions [War 07] . 

Directed last passage percolation. "Directed" refers to the constraint of 
not being allowed to turn back in time. In the physics literature one speaks 
of a directed polymer in a random medium. Shape theorems are proved, 
e.g., in |Kes86] . While the issue of fluctuations had been repeatedly raised, 
sharp results had to wait for [JohOOj and |PS02bj . Growths models naturally 
lead to either point-to-point, point-to-line, and point-to-random-line last 
passage percolation. For certain models one has to further impose boundary 
conditions and/or symmetry conditions for the allowed domain. In (IMl) one 
takes the max, thus zero temperature. There are interesting results for the 
finite temperature version, where the energy appear in the exponential, as 
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in | CY06] . 



KPZ equation. The seminal paper is |KPZ86] . which generated a large 
body of theoretical work. An introductory exposition is |BS95] . The KPZ 
equation is a stochastic field theory with broken time reversal invariance, 
hence a great theoretical challenge, see, e.g., |Las98j . 

Review articles. A beautiful review is |Joh06] . Growth models, of the type 
discussed here, are explained in |FP06|lFfer08] . A fine introduction to random 
matrix techniques is |Sas07j . |KK08] provide an introductory exposition to 
the shape fiuctuation proof of Johansson. The method of line ensembles is 
reviewed in |Spo06| . 
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